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Abstract 

Eugenio Bianchi and Matteo Smerlak have found a relationship between 
the Hawking radiation energy and von Neumann entropy in a conformal field 
emitted by a semiclassical two-dimensional black hole. We compare this rela¬ 
tionship with what might be expected for unitary evolution of a quantum black 
hole in four and higher dimensions. If one neglects the expected increase in the 
radiation entropy over the decrease in the black hole Bekenstein-Hawking A/4 
entropy that arises from the scattering of the radiation by the barrier near the 
black hole, the relation works very well, except near the peak of the radiation 
von Neumann entropy and near the final evaporation. These discrepancies are 
calculated and discussed as tiny differences between a semiclassical treatment 
and a quantum gravity treatment. 
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1 Introduction 


Eugenio Bianchi and Matteo Smerlak mq have found a beautiful formula relat¬ 
ing the Hawking radiation energy flux F and the retarded time derivatives of the 
von Neumann entanglement entropy S at future null inhnity for a two-dimensional 
conformal held theory in a hxed two-dimensional classical or semiclassical spacetime: 

Here F is the energy hux at future null inhnity (X+) as a function of the retarded 
time M, S is the renormalized entanglement entropy of the radiation at up to 
the time u, c is a constant that depends on the conformal held, and an overdot 
represents a derivative with respect to the time u. This formula has been applied to 
many solvable models of gravitational collapse by Bianchi, De Lorenzo, and Smerlak 
0 . 

Here we wish to compare the predictions of this formula with what is expected to 
be the case for the energy hux and the von Neumann entropy of Hawking radiation 
of massless helds from a four-dimensional spherically symmetric black hole. For the 
emission of a conformally invariant scalar held, one might expect the Hawking radi¬ 
ation to be dominated by the scalar held modes that have zero angular momentum 
(S-waves) and hence are spherically symmetric, ehectively reducing the problem to a 
two-dimensional one for which one might expect the formula of Bianchi and Smerlak 
to apply, at least to some level of approximation. 

We examine whether this is indeed the case. We hnd that during most of the 
Hawking emission by a black hole of initial mass Mq that is large in Planck units, 
the hrst term on the right hand side of Eq. Q dominates over the second term, and 
the power-law dependence of this hrst term on the mass M is the same as that of the 
hux F on the left hand side. Therefore, by choosing the constant c appropriately, 
one can get a good match between the left and right hand sides. 

However, there are some caveats for this description. 

First, for a black hole that starts in a pure state and evaporates away completely 
by a unitary process that does not lose information, the von Neumann entropy S 
of the radiation starts at zero at the beginning of the evaporation (when there is 
not yet any radiation, so that its entropy is zero) and goes back to zero at the end 
of the evaporation (when all of the information is in the radiation, so that it is in 
a pure state again with zero von Neumann entropy). In between, the radiation is 
entangled with the black hole, so that its von Neumann entropy S{u) rises from zero 
to a maximum (at what has been called the ’’Page time” and which we shall denote 
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by M = Mp, at which retarded time we shall say the entropy is S{up) = Sp and the 
black hole horizon radius is R{up) = Rp) and then drops back down to zero. 

Because the scattering of modes by the potential barrier near the black hole 
horizon makes the evaporation a nonadiabatic process, during the stage in which 
S{u) is rising, its increase is greater, by a factor h > 1 that is a constant for the 
massless helds that we shall consider here, than the decrease in the Bekenstein- 
Hawking entropy S'bh = ^/(4G), where A is the area of the black hole event horizon. 
(For example, for a four-dimensional Schwarzschild black hole emitting into empty 
space, numerical calculations [1] have shown that for the emission of just massless 
photons and gravitons, b ~ 1.4847.) The Bekenstein-Hawking entropy is believed to 
give the leading term for the logarithm of the number of quantum states of a black 
hole of area up to A and hence a good approximation for the maximum von Neumann 
entropy a black hole of area A can have. Once the effective number of accessible 
quantum states for the Hawking radiation exceeds that of the black hole (after Up), 
the von Neumann entropy of the radiation is limited by S'bh and is expected to be 
very close to this limit 0 E]. Therefore, before Up, one expects that S ~ —bSBu, 
whereas after Up , one expects that S ~ +Sbh- On the other hand, the energy flux 
F is proportional to *Fgjj, with the same coefficient, for both time periods. Then for 
Eq. ([^ to be valid for both periods, c would have to decrease by a factor of 1/6^ 
when the retarded time u crossed Up, which is not consistent with the assumption 
that c is a constant depending only on the massless helds being emitted. 

Second, when S{u) goes from increasing before Up to decreasing after Up, there is 
a brief period of retarded time u during which the second term (with the second time 
derivative of the von Neumann entropy, S) on the right hand side of Eq. ([^, which 
is negative, dominates over the positive semi-dehnite hrst term (with the square of 
the hrst time derivative of the von Neumann entropy, 3“^), leading to a negative 
expression for the energy hux F. If one trusted Eq. ([^ during this stage of the 
evaporation, it would seem that the black hole would have to gain a bit of energy 
at this time, rather than continuously losing energy. 

Third, if one extrapolates the semiclassical approximation for 3 ~ +*Fbh as a 
function of the black hole mass M down to Planck and sub-Planck values (where 
the semiclassical approximation is not believed to be valid), one again gets a regime 
in which the second term on the right hand side of Eq. ([^ dominates over the hrst 
term and again gives a negative expression for the energy hux. 

Here we shall discuss these caveats that challenge the validity of Eq. ([^ for a 
quantum four-dimensional black hole. The hrst caveat is perhaps the most serious, 
but if it can be swept under the rug, then the second caveat just gives a very small 
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violation of monotonic mass loss and might be an artifact of the approximation 
in Eq. Q of a single dehnite (semiclassical) metric. The third caveat of another 
violation of monotonic mass loss near the end of the Hawking evaporation can be 
avoided even more simply by a slight modihcation of the time dependence of the 
von Neumann entropy near the hnal decay of the then-tiny black hole. 

One might object (as an anonymous referee has) that there is a priori no rea¬ 
son to expect two-dimensional gravity models to match higher-dimensional gravity 
quantitatively, so it should not be surprising when they do not. 

Indeed (to go beyond the brief objection of the referee), one might give a simple 
reason why something like Eq. ([^ holds approximately in d = 2 spacetime dimen¬ 
sions but not in higher dimensions. The hrst term on the right hand side of Eq. ([^ 
has an extra power of the entropy S than the second term, whereas both terms have 
the same number of time derivatives (two). Therefore, when the entropy is large and 
when the square of the hrst time derivative of the logarithm of the entropy is not 
much smaller than the second time derivative of this logarithm, one would expect 
the hrst term to dominate, and then Eq. ([^ says that the ratio between the energy 
hux F (the power or luminosity) and the square of the entropy hux S should be a 
constant independent of the gravitational system. 

For objects of temperature T and hxed shape with linear size R that is much 
larger than a thermal wavelength d spacetime dimensions, massless radiation gives 
an energy hux that goes as F ~ whereas the entropy hux goes as ^ ~ 

Therefore, S'^/F ~ {RTY~‘^. In d = 2 the right hand side is a constant 
independent of the size R and the temperature T of the object, so one would indeed 
expect that the energy hux would be roughly proportional to the square of the 
entropy hux with a constant of proportionality depending on the number and types 
of helds in the thermal radiation. 

Now for black holes of hxed shape (e.g., spherical, or else rotating with one or 
more hxed dimensionless rotation parameter that determine the shape), the only 
quantity that sets the scale is the size R, which determines the temperature T to be 
some size-independent constant divided by R, so for asymptotically hat black holes 
RT is just a shape-dependent constant. Therefore, in this special case S'^/F ~ 
{RTY~^ for massless radiation is indeed independent of the size of the black hole. 

However, if another relevant scale is present, then RT is not necessarily inde¬ 
pendent of the size of the black hole, and so S'^/F can depend on the size. For 
example, for spherical black holes in asymptotically anti-de Sitter spacetime that 
has an additional scale i given by the cosmological constant, RT depends on R/i, 
and for i? S> £, S'^/F ~ {R/iy~^ [7j, which is not independent of the hole size R. 
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Despite these objections to applying the two-dimensional Bianchi-Smerlak Eq. 
([^ to higher-dimensional black holes, it may be of interest to examine the one ideal¬ 
ized case in which it might be hoped to apply, spherical black holes in asymptotically 
flat spacetime so that the only length scale is that of the black hole and RT ~ 1, 
independent of the one size parameter for such black holes. 


2 Hawking Radiation Energy and Entropy 


We shall generally use units in which the speed of light c, the reduced Planck con¬ 
stant % = h/ (27r), and the Boltzmann constant ks are set to unity, c = h = ks = 
but not the Newtonian constant of gravitation G = 6.67384(80) x 10“^^m^kg“^s“^ = 
2.61210(31) X 10-™m2 = 2.90635(35) x 10-®V = 2.11095(25) x lO^^kg"^. How¬ 
ever, sometimes we shall refer to Planck units, in which one sets G = 1 as well as 
c = h = ks = I- 

Let us consider a d-dimensional, non-rotating, spherical, uncharged Schwarzschild- 
Tangherlini black hole of horizon radius R, with the metric 

= -(1 - + (1 - + ^^d^fd-2), ( 2 ) 

where /i = R‘^~^ and is the metric of a (d — 2)-dimensional unit sphere. The 

Arnowitt-Deser-Misner (ADM) mass M of the black hole is given by 


M 


(d — 2)D(rf_2)p 

IQtiG 


(d - 2)D(,_2)i?"-^ 
IQnG 


(3) 


where 


^{d 2) (4) 

is the area of a unit (d — 2)-dimensional sphere. The area of the black hole horizon 
is 


A — il(d-2)R'^ ^5 

and the Bekenstein-Hawking thermodynamic entropy is 

A dvr / WnG \ ^ ^ 

4G 4G d-2 V(d-2)D(rf_2); 


(5) 

( 6 ) 


Introducing the null retarded/outgoing coordinate u 
vanced/ingoing null coordinate v = {t + r*), where 


dr* 


dr 




r 

rpd—Z 


•> 


{t — r*) and the ad- 


(7) 
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we can rewrite the metric ([^ in the following form: 


ds'^ 



v) J 


dudv + r^dujf^_ 2 y 


( 8 ) 


Let us consider now that this black hole is evaporating by the process of Hawking 
radiation to purely massless radiation. Then the black hole mass M, horizon radius 
parameter /i = horizon area A = and Hawking-Bekenstein 

entropy ^bh = ^/(4G') are all functions of the retarded time u. For a spherical 
uncharged black hole large in Planck units {GM‘^ S> 1), the Hawking radiation 
energy flux of the black hole is 


dM a a / {d — 

du B? /i2/('^-3) IbvrGM / ’ 


where a is a dimensionless constant which can be found from the Hawking emis¬ 
sion rate of massless particles. Numerical calculations [HI El EHl [m E] show that 
for a four-dimensional black hole, the emission power is approximately 3.3538 x 
in photons and 0.3836 x 10“^/(GM)^ in gravitons, for a total emis¬ 
sion of massless radiation energy about 3.7474 x This gives a ~ 

0.00014990 ~ 1.5 X lO'L 

A black hole of initial mass horizon radius Rq and initial mass Mq = 

(d — 2)r2(d_2)i?o“^/(167rG) (say at m = 0), emitting purely massless radiation and 
hence having a Hawking energy flux F = a/R^, will then last a lifetime 


Ut = 


(d 2)(d 3)r2(rf_2) _ d 3 


pa—i _ 

-flO — 


IGttG 


2 

d-Z 


167r(d — l)Ga {d — l)a \{d — 2)^l(^d-2), 

Then the black hole mass M{u) as a function of u is given by 


d-1 


( 10 ) 


Miu) = Mo 1- 

Ut 


u 


d-Z 

d-1 


( 11 ) 


The semiclassical approximation for the black hole emission gives the time-dependence 
of the coarse-grained entropy of the black hole and thus of the emitted radiation as 


*S'rad(M) — &[5 'bh(0) — 5'bh(m)] 


&*S'bh(0) 



( 12 ) 


where S'bh(O) = Ao/(4G) = r2(d_2)i?o“^/(4G) is the initial Hawking-Bekenstein 
entropy of the black hole with mass Mq, S'bh('w) is the Hawking-Bekenstein entropy 
of the black hole when its mass is M{u), and b is the ratio by which the increase in 
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the coarse-grained entropy of the Hawking radiation (e.g., ignoring its entanglement 
with the black hole) is greater than the decrease in the coarse-grained Bekenstein- 
Hawking entropy of the black hole. Numerical calculations mi show that for a 
massive 4-dimensional Schwarzschild black hole that emits essentially only photons 
and gravitons, b ^ 1.48472. 

We assume that the semiclassical Bekenstein-Hawking entropy of a nonrotat¬ 
ing uncharged black hole is a good approximation for the maximum von Neumann 
entropy of a black hole of the same energy, at least if one neglects the entropy as¬ 
sociated with the location and/or motion of the black hole in a space sufficiently 
large that this could in principle rival the Bekenstein- Hawking entropy. Second, 
the semiclassical entropy calculated for the Hawking radiation is the maximum von 
Neumann entropy for radiation with the same expectation value of the number of 
particles in each of the modes. Under the extra assumptions that the black hole 
starts in essentially a pure state, that the Hawking evaporation is a unitary process, 
and that we can neglect the interaction with other systems, the von Neumann en¬ 
tropy of the evaporating black hole equals that of the Hawking radiation that has 
been emitted. 

The Bekenstein entropy of the black hole decreases monotonically with time and 


the radiation entropy (12) increases monotonically with time. The two values cross 


at the so-called ‘Page time,’ at 


Up — Uf 


1 - 


d-1 

b \ <1-2 


b + 1 


(13) 


Then when one applies the assumption of unitarity to the results of the semiclas¬ 
sical approximation, the von Neumann entropy of the radiation during the entire 
evaporation is expected to be given to a good approximation by the minimum of 
the dimensions of the radiation and black hole subspaces of the Hilbert space |1] , so 


s ^ 5'approx = &[5'bh(0) - Sbr{u)]9{up - u) + Sbh{u)9{u - Up 


— &*S'bh(0) 


1 - — 
Ut 


d-2- 

d-1 


U 


- m)S'bh(O) ( 1 - — ) 9{u-Up), (14) 


d-2 

d-1 


where 9 is the Heaviside step function. The maximum of the von Neumann entropy 
is at M = Up, where the coefficients of the two step functions above coincide. 
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3 Comparison with the Bianchi-Smerlak Formula 

Now let us compare the approximations above for the energy flux, Eq. @ , and 
entropy, Eq. (12), of the Hawking radiation with the Bianchi-Smerlak formula ([^. 

The Hawking temperature for the d-dimensional spherical uncharged Schwarzschild- 
Tangherlini black hole of horizon radius R is the surface gravity k divided by 2tt. 
With —goo = = V = 1 — = 1 — and with the radial proper 

length element dl = ^/^dr, the surface gravity is the following expression evaluated 
on the horizon, r = R: 


K = 


d^/^0 d{V^/^) IdV d-3R 


d-3 


This then gives 


dl V-R^dr 2 dr 


^ K d — 3 

T = — =- 

27r 4:7tR 


„d—2 


(15) 


(16) 


Now the hrst law of black hole thermodynamics gives dM = RdS^n, so M = 
—F = —a/R^ gives 


A _M _ . _ 

*->BH — — — —— — 


dvra 


T 


T 


{d-3)R' 


From M = {d — 2)VL(^^_2)R^ ^/(IbvrG), one also gets 


and hence 


R = 


‘S'bh — 


TOttGo 


(d-2)(d-3)H(,_2)i?'='-2 

~(d-2){d-3yn(^d-2)R'^' 


Then from Eq. (14), one gets 


S^S: 


approx 




Anab 


{d-3)R 


6{up — u) — 


Ana 


{d-3)R 


9{u - Up), 


and 




QAn^Ga^b 


approx (d-2){d-3)m^d-2)R^ ^ 

GAn^Ga^ 


6{up — u) 


{d-2){d-3)^^l(d-2)R^ 


at ^ 47ra(6 + l) 

dya — Up) ——— d[u — Up). 


{d — 3)R 


(17) 


(18) 


(19) 


( 20 ) 


(21) 


V 

















Plugging this into the Bianchi-Smerlak formula ([^ then gives 


F^Fbs = -i-S^ 


271 \c ‘S'approx 


487ia%‘^ 


+ 


32TiGa^b 


c{d-?>fR^ {d - 2){d - ?>fni^d-2)R^ 

487ra^ 827iGa? 


c{d-?>fR^ {d - 2){d - ?>fn(d-2)R‘^ 

2a(6 + 1) 


6{up — u) 

6{u — Up) 


{d - 3)Rj 


-6{u — Up). 


( 22 ) 


For u < Up, the ratio of the second term (proportional to S) to the hrst term 
(proportional to S"^) is c/[66((i —2 )S'bh]; for u > Up it is —c/[6((i —2 )S'bh]- Therefore, 
except very near u = Up (the so-called Page time when the von Neumann entropy 


of the radiation reaches its maximum), where Eq. (14) is not a good approximation 


to the von Neumann entropy S that actually rounds off rather than having a kink 
(see below), the second term is much smaller than the hrst term whenever the black 
hole is large in Planck units, so that Fbh 3> 1, which breaks down near u = Ut, near 
the hnal evaporation of the black hole. That is, except for the retarded time u near 
either Up or ut, 


F Fbs 


487ra^ 


c(d-3)2i?2 


\jj^6{up — u) + 6{u — Up)]. 


(23) 


This would ht well with what one expects for a large black hole emitting almost 
entirely massless radiation, F = a/R^ hy Eq. (§ if 


c = 


48™ 

(d-3)‘ 


\b‘^6{up — u) + 6{u — Up)]. 


(24) 


However, c is supposed to be a constant, not changing between the periods u < Up 


and u > Up, which is in conhict with Eq. (24) if 6 7 ^ 1. Therefore, the Bianchi- 


Smerlak formula ([^ does not seem to work well for a black hole evaporating so that 
during the early stage the von Neumann entropy of the radiation is larger by a factor 
6 > 1 than the decrease in the Bekenstein-Hawking entropy S'bh = 24/{4G) of the 
black hole (not the von Neumann entropy during this stage, though it is to a good 
approximation for u > Up). This does seem to be a serious problem with applying 
the Bianchi-Smerlak formula (actually derived just for emission in two dimensions 
from a hxed spacetime metric) to black holes in d > 3 dimensions with a quantum 
spacetime. 
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See Figure 1 for a plot of M{u), the time integral of —F = —a/R'^ and of —-Fbs 
for the two constant values of c, neither of which matches that of the expected 
-F = -a/R?. 



Vjf 


Figure 1: Plot of black hole mass vs. time in 4-dimensions. Solid black line is 
the expected M{u) for dM/du = —F = —a/(2M)^. Dashed blue line is M{u) for 
dM/du = —Fbs with c = 487ra. Dashed-dotted red line is M{u) for dM/du = —Frs 
with c = 487ra6^. 


Nevertheless, let us proceed by sweeping this problem under the rug and pro¬ 
ceeding to the next issue, which is how the Bianchi-Smerlak formula for the flux 
behaves near the peak in S'approx, which is at u = Up. If we use Eq. (14) for the 
entropy S, this has a kink in it (even if we set 6 = 1) at u = Up, so the resulting 
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approximation (22) for the flux F has a negative delta function at u = Up, giving a 
sudden increase in the black hole mass. In reality S{u) will be smooth near u = Up, 
so let us look for an improved approximation giving such a smooth S{u). 

In O E], it was noted that if one has a random pure state (using the Haar 
measure) in a Hilbert space that is the tensor product of two subsystem Hilbert 
spaces of dimensions m and n, both large, then for m <n the average von Neumann 
entropy of each subsystem is 

Tfl 

S Inm-. (25) 

2n ^ ’ 


For a black hole that is initially in a pure state surrounded by vacuum (another 
pure state), we can take m and n to be the effective Hilbert space dimensions 
of the black hole (using the exponential of the Bekenstein-Hawking entropy S'bh = 
A/{AG) as an estimate of the dimension of the Hilbert space of black holes of horizon 
areas close to A) and of the radiation (using the exponential of the thermodynamic 
radiation entropy Srad = ^[*S'bh(0) —S'bh(w)] as an estimate of the effective dimension 
of the Hilbert space of the radiation that would be entangled with the hole; if the 
black hole is emitting into asymptotically flat empty space, the total Hilbert space 
dimension for the helds in this space would be inhnite, even for bounded total energy, 
but nearly all of that Hilbert space would develop negligible entanglement with the 
black hole and hence can be assumed to be in nearly a pure vacuum state and not 
contribute much to the entanglement entropy). 

For brevity in the equations below, let the Bekenstein-Hawking thermodynamic 
entropy of the black hole be 


s = Sbb{u); So = 5'bh(0). 


(26) 


Then the radiation entropy is approximately b{so — s). 

Therefore, for u < Up, the effective dimension of the smaller subsystem (the 
radiation) is m ^ exp {bso — bs), and that of the larger subsystem (the black hole) 


is n exp (s), whereas for u > Up, m and n are reversed. Eq. (25) then gives 


^ ~ *^smootn 


bso — bs — - exp {bso — bs — s) 


+ 


s — - exp (s -|- 6s — 6so) 


9{up — u) 

6{u — Up). 


It is convenient to define 


(27) 


X = bso - {b + l)s - (^; + l)(sp - s), 


(28) 
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where 


Sp 


6 + 1 


So 


(29) 


is the value of Inm = Inn where they are equal, at u = Up. This gives x < 0 for 
u < Up, and x > 0 for u > Up. Then Eq. (27) gives the von Neumann entropy as 

b 


S 


Sp + -—-X - -e" 6{-x) + s. 


’p 6 + 1^ 2^ 


' b + r 2 

The hrst two time {u) derivatives of the von Neumann entropy are 

6 1 


6{x) 


where 


5 


+ 




1 1 -X 

-1—e 

6 + 1^2 




6 + 1 


2^ I ^ 


9{-x) + 


+ 


1 1 -X 

-1—e 

6 + 1^2 


X 


i = _(t + 1), ^ _(6 + D^SMu) = -"jz-sjr 

i = _(t + 1)S ^ _(6 + 1 )Tsb„(«) = 


(30) 

(31) 

0(x), 

(32) 

(33) 

(34) 


One can see that S, S, and S are all continuous at u = Up, but the third derivative, 
d^S/du^, is not. Also, one can see that, except near the hnal stage of the Hawk¬ 
ing radiation, when the area A gets near the Planck value that is the Newtonian 
gravitational constant G in our natural units with h = c = ks = I, the Bekenstein- 
Hawking entropy s is very large, so |x| <C x^, meaning that we can effectively ignore 
the X terms. 

Plugging these equations into the Bianchi-Smerlak Eq. ([^ then give the flux 

2 


26 


dvrc \ 6 + 1 


1 x^ 


^e"x2 + ^ 


dvr 


26 


dvrc \ 6 + 1 


X 


1 -X.2 , 1 

—e X + 


dvr \ 6 + 1 
2 


X 


dTT 


dTT 


6 + 1 


+ e" 


X 


9{—x) 

9{x). (35) 


Because of the continuity of S and S, the flux F = —M = —dM/du is also contin¬ 
uous at M = Up, but F = dF/du = —d'^M/du^ is not. If we do drop the x terms in 
the expression above for F, we get 




x^ 

dvr 


3/26 
c V6 + 1 



— e 


9{-x) + 


3 / 2 
c + 1 
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See Figure 2 for a plot of the normalized Bianchi-Smerlak flux {R^/a)FBs, for 
d = A and c = dSvra, versus S'bh('Wp) — Sbh{u) = x/{b + 1). 


—F- 

a a du 



Sbh{up) — Sbh{u) 


Figure 2: Plot of Bianchi-Smerlak flux formula near the peak of S. Solid line is 
(i?^/a)FBs for d = A and c = dSvra. Dashed red line is the asymptotic value for 
u — Up —R, a factor of 6^ too large with this constant value of c. Dotted blue 
line is the asymptotic value for u — Up ^ R, matching the expected value of the 
normalized flux. 
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Since S{u) is turning over from increasing to decreasing, there is always a range 
in X = bsQ — {b + l)s for which the {6/c)S‘^ term in the Bianchi-Smerlak Eq. Q for 
the flux F is smaller in magnitude than the negative S term, so F is negative in this 
range. For [{b — l)/(6 + 1)]^ > c/3, which is the case for large 4-dimensional black 
holes emitting just photons and gravitons, the range is entirely in x > 0 and is 


In 



2 In (a/1 -|- 5 -|- \/6) < X < In 



2 In (a/1 (5 v^). 


(37) 


where 


5 = 


(6 -|- l)c 27ia{b -|- 1) 
24 ^ (d-3)2 


0.0023402, 


(38) 


using Eq. (24) for c = 487ra/((i — 3)^ for u > Up where the x > 0 region with 
negative F occurs, and then at the end (as shall also be done in successive equations 
below) inserting the numerical values above for a and for b for a large 4-dimensional 
Schwarzschild black hole that emits only photons and gravitons. 

If we dehne 


r/= 41n (vT+~5-|-a/ 5) = 4sinh ^ a/ 5 Ri 0.19343, (39) 

then integrating the negative flux over the time period [{d — 3)yRp\/[4:7ia{b + 1)] ~ 
AlRp over which it occurs (with Rp being the radius of the black hole horizon then) 
gives a mass increase during this period of 

AM = Tp ~ 0.0004863Tp, (40) 

where Tp = (d —3)/(47ri?p) is the Hawking temperature of the black hole at the time 
very near u = Up when the von Neumann entropy S is maximized and when the flux 
F = —dM/du is negative, leading to the increase in the mass M. 

We see that according to the Bianchi-Smerlak formula Eq. ([^, the mass would 
indeed increase, but the increase would be extremely small for a large-mass black 
hole in 4 dimensions, less than 0.05% of the energy Tp of one thermal photon at the 
Hawking temperature. 

Furthermore, because of random root-A fluctuations in the number of particles 
emitted, which in 4-dimensions is y/N ~ R/y/G, with each particle taking a time 
~ i? to be emitted, one might consider the Hawking evaporation to give an ensemble 
of semiclassical spacetimes, each with its own time at which the von Neumann 
entropy of the radiation is maximized. The rms spread in these times would be 
of the order of ~ R^/yfG, which for a solar-mass black hole of about 10^® Planck 
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masses would be of the order of 10^® Planck times or 10^® times the present age 
of the universe. If the Bianchi-Smerlak formula applied to each spacetime in the 
ensemble, then there would be a negative mass flux of energy ~ —0.0005 Tp over a 
time less than a millisecond but uncertain by a time of the order of 10^® times the 
duration of the negative flux, a time roughly hfteen orders of magnitude larger than 
the present age of the universe. Therefore, the effect of this negative energy flux, if it 
really occurs, seems to be quite negligible and virtually impossible to detect. Surely 
it would be washed out by any averaging over the uncertain time at which a random 
member of the ensemble of radiating black hole spacetimes reached the point where 
one might expect the von Neumann entropy of the radiation to be maximized. 

One can also ask how small one can make the negative energy emitted in a two- 
dimensional black hole mirror model in which one has control of the motion of the 
mirror and hence of the time evolution of the entropy, S{u). In the four-dimensional 
black hole example above in which S{u) was the specihc function given by Eq. (30), 
we found that the negative energy emitted was just a small fraction (but still of 
the general order of magnitude of unity) of the energy of a single thermal quantum. 
However, by choosing the mirror trajectory appropriately for a two-dimensional 
model, we can get the negative energy emitted to be of the order of a thermal 
quantum divided by the maximum von Neumann entropy and hence enormously 
smaller than one thermal quantum. For example, consider the following model with 
four retarded time periods, each of length r: 

(1) For —2r < u < —r, let S{u) = Si{u/t + 2)^, which goes from 0 at u = —2r 
to Si ai u = —T. S = {2Si/t){u/t + 2) goes from 0 at u = —2r to 2Si/t at m = —r. 
S = 2S'i/r^ is a constant during this period. 

(2) For —r < n < 0, let S{u) = 51-1- (c/6) In [1-|-(125i/c)(M/r-|-1)], which 

goes from Si at u = —r to Si + (c/6) In (1 -|- 125i/c) at u = 0. S = (25i/r)/[l -f 
(125i/c)(M/r -|- 1)] goes from 2Si/t at u = —r to (25i/r)/(l -|- 125i/c) at u = 0—. 
S = -[24(5i/r)Vc]/[l (125i/c)(M/r 1)] = -(6/c)5^. 


Now one makes the entropy a symmetric function of the retarded time u, so 
S{u) = S{—u). At u = 0, S jumps from (25i/r)/(l -|- 125i/c) at u = 0— to 
— (25i/r)/(l -|- 125i/c) at u = 0-)-, so to cover the period —r < u < t during 
which 5^ is continuous, S = —(6/c)5^ — (45i/r)/(l -|- 12Si/c)5(u). At the other 
two junctions between the four regions, both S and S are continuous functions of 
M, though S is not. 

(3) For 0 < M < r, S(u) = 5i -|- (c/6) In [1 — (125i/c)(M/r — 1)]. 

(4) For r < n < 2r, S{u) = Si{u/t — 2)^. 
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Now the Bianchi-Smerlak Eq. Q gives the energy flux 


as 


F = 


12Sf f\u\-2T 


ncT^ 


+ 


Si 


vrr^ 


9(2r-H)9(H-r)-^-i^^%), (41) 


which is positive for r < |m| < 2r, zero for 0 < |m| < r, and has a negative 
delta function at m = 0. Positive energy Ei + E 4 = { 8 Sf + 2cSi)/( tict) is emitted 
during the periods (1) and (4) with t < \u\ < 2r, and negative energy —^negative = 
—4cS'i/[27rr(12S'i + c)] is emitted at m = 0 when S jumps downward. 

When Si is very large in comparison with 1, c, and c^, the maximum von Neu¬ 
mann entropy is Smax = >S'i-|-(c/6) In (1 -|- 12S'i/c) ~ Si, and the total energy emitted 
is E = EiEE/1—E 

negative (8S'^-|-2cS'i)/(7rcr)—4cS'i/[27rr(12S'i-|-c)] ~ 8S'^/(7rcr) ~ 
8>5 '^ay /( ttct), which is approximately a factor of 16/3 larger than the minimum en¬ 
ergy E min = 65 'ma^ / (ttcAm) that can be emitted during a period of retarded time 
Au (here 4r) if S and S both start at zero at the beginning of the period and end 
at zero at the end of the period. (The minimum is attained by having S increase 
uniformly from zero to S'max during the first half of the total retarded time period 
Au and then decrease uniformly back to zero during the second half of the period.) 

However, the magnitude of the negative energy emitted in the example above 
is only ^negative c/(67rr), with no factor of S'max, so it is smaller by a factor 
of the order of l/^max than the magnitude of the positive energy emitted. To 
put it another way, one might say that an average quantum emitted has energy 
.E/S'max ~ 8S'max/(7rcr), and in comparison, Enegative/ci; ~ c^/(48S'i) < 1. That 


u 


is, in the model above, the negative energy emitted is much smaller than the average 
energy of a quantum emitted, smaller by a factor of very roughly the reciprocal of 
the maximum von Neumann entropy. In this sense the negative energy flux implied 
by the Bianchi-Smerlak Eq. ([^ when the von Neumann entropy S{u) is at or near 
its maximum value can be made quite negligible if this von Neumann entropy is 
very large. 

Now let us go back to the four-dimensional black hole model and turn to the 
final stages of its Hawking evaporation. If one takes the approximate formula Eq. 


(14) for the von Neumann entropy that is valid except near u = Up and inserts it 
into the Bianchi-Smerlak formula Eq. ([^ for the flux to get Eq. (22), one finds that 
one also gets a negative flux near the end at the black hole lifetime, at u ~ Uf. 


However, one would not expect that Eq. (14) would be valid when the black hole 


has evaporated down to near the Planck size. One does not know what the actual 
expressions should be (which also should take into account the uncertainty in the 
time Ut, analogous to the uncertainty in the time Up, that would be of the order of 
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Ut/\/ *S'bh( 0) V^, an expression valid in any number of spacetime dimensions 

greater than 3 where we can have black hole evaporation). Nevertheless, one can 
easily come up with expressions that give a positive energy flux all the way down to 
zero mass and size for the black hole. 

For simplicity, we shall just consider d = 4. If one dehnes X = R‘^/{8Ga), then 
for a large Schwarzschild black hole, which has X S> 1, the Bekenstein-Hawking 
entropy is Reh = = SnaX, and the Hawking evaporation gives 


R = 2GM = 


2Ga 


1 

'ax' 


(42) 


If now one keeps this Eq. (42) for R but replaces the Bekenstein-Hawking entropy 
formula by 


5 = 


8naX 


C2 


(43) 


one 

flux of the form 


F = — 


1 -l- 1/X Rbh 4“ Svra 

can calculate that the Bianchi-Smerlak formula Eq. ([^ with c = 487ra gives a 

_ + 8V^ + 9X 1 

2{X + !)■> J ' ’ 

This never goes negative but instead increases monotonically with the time u as 
R decreases, always remaining within a factor of 3 of the expression F = a/R? 
from Eq. Therefore, there is no problem constructing S{u) that is close to the 
Bekenstein-Hawking black hole entropy Feh = ^/(4G) for a large hole but deviates 
from it when the hole gets small so that the Bianchi-Smerlak formula for the flux 
stays positive no matter how small the black hole gets. The flux can even have the 
asymptotic form F = a/R^ for very small R? jG = 8aX as well as for very large 
R? jG, as the example above shows. 


See Figure 3 for a plot of the factor in the square brackets in Eq. (44), {R^/a)F, 
versus the retarded time u — Ut, near the end of the black hole evaporation. 
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Figure 3: Plot of the normalized Bianchi-Smerlak flux near the end of the Hawking 
radiation: {R‘^/a)F versus retarded time u — Ut, for S = S'|j£/(S'bh + Svra) instead 
of for S = S'bh- 
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4 Conclusions 


In conclusion, the Bianchi-Smerlak formula Eq. Q, which was derived for a 2- 
dimensional model black hole with quantum held theory on a dehnite spacetime 
metric (say with boundary conditions given by a moving mirror), seems to work fairly 
well for quantum black holes in higher dimensions, but there are some inadequacies 
in the model. Perhaps the most serious is the fact that one cannot match the entropy 
and energy hux both before and after the peak in the von Neumann entropy with 
a hxed value of the constant c in the Bianchi-Smerlak formula if the black hole 
evaporation process is not adiabatic but generates extra entropy in the radiation in 
excess of the decrease in the Bekenstein-Hawking entropy A/ (4G) of the black hole. 

Another potential problem is the negative flux predicted near the peak value 
of the von Neumann entropy S', which would cause the black hole briefly to gain 
mass instead of losing it monotonically. However, using the numerical results for 
Hawking emission of the known massless particles in 4-dimensions leads to a very 
tiny increase in the mass, less than 0.05% of the energy of a single quantum of the 
energy of the Hawking temperature of the black hole at that time. One can further 
argue that this tiny effect would be totally swamped by the quantum uncertainty 
of when in each semiclassical spacetime the von Neumann entropy would reach its 
peak. 

A third potential problem is the fact that extrapolating the expected behavior 
of the von Neumann entropy of a large black hole down into the hnal stages in 
which the black hole gets near the Planck mass leads also to a negative flux of 
energy. However, unlike the negative flux at near the peak in S', which cannot be 
avoided in the Bianchi-Smerlak formula, one can modify the assumed behavior of 
the von Neumann entropy near the hnal evaporation so that the Bianchi-Smerlak 
formula continues to give a positive and monotonically rising hux at the end of the 
evaporation. 
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